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Abstract: We study the interaction of particles with a domain wall at 
a symmetry-breaking phase transition by perturbing about the domain wall 
solution. We find the particulate excitations appropriate near the domain 
wall and relate them to the particles present far from the wall in the uni- 
form broken and unbroken phases. For a quartic Higgs potential we find 
analytic solutions to the equations of motion and derive refiection and trans- 
mission coefficients. We discover several bound states for particles near the 
wall. Finally, we apply our results to the electroweak phase transition in the 
standard model. 
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1 Introduction 



During a first-order phase transition in which a gauge symmetry is spon- 
taneously broken, a number of effects depend on the scattering of particles 
from the domain wall between the phases of broken and unbroken symmetry. 
For instance, at the electroweak phase transition, a difference in the reflec- 
tion coefficients of particles carrying opposite quantum numbers could play a 
role in producing the baryon asymmetry of the universe. Analytic solutions 
for quark scattering from domain walls of two different profiles are already in 
the literature , but so far there has been no analogous treatment for the 
bosons of the theory. The study of boson scattering at a symmetry-breaking 
phase transition is complicated by the fact that the two phases have differ- 
ent particle contents, as explained by the Higgs mechanism. Addressing this 
difficulty is the main focus of the present work. 

We begin by considering the simplest gauge theory with spontaneous 
symmetry breaking, the Abelian Higgs model. To permit a static domain 
wall to exist, we study the theory at its transition temperature. We find the 
modes of particulate excitation appropriate near the domain wall and express 
the interaction of these modes with the wall in terms of a one-dimensional 
scattering potential. We relate these internal modes to the particle modes 
present in the asymptotic broken and unbroken phases. All these results are 
valid for a general Higgs potential. 

In section ^ we specialize to a quartic Higgs potential and obtain ana- 
lytic solutions to the one-dimensional scattering equations for the internal 
modes, including solutions which describe bound states. Using the connec- 
tions between internal and asymptotic modes derived in section we compute 
scattering probabilities for the asymptotic particle states. 

Having completed our analysis of the Abelian Higgs model, we proceed 
in section ^ to study the electroweak phase transition in the standard model. 
As a nonabelian gauge theory with partial symmetry breaking, the standard 
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model should be representative of the entire class of spontaneously broken 
gauge theories. We find that our results for the Abelian Higgs model may be 
adapted without difficulty to the case of the bosons in the standard model. 
We complete our discussion of the standard model by applying the methods 
of the previous sections to the case of the fermions. Specifically, we obtain in- 
ternal and asymptotic modes, scattering potentials, and connection matrices 
for fermion scattering. 

2 Internal and asymptotic modes 

2.1 The model 

The Abelian Higgs model contains two fields, a complex scalar field and 
a C/(l) gauge field Ai^. The Lagrangian is 

L = -\f,,F^'' + {D^<t^r{D^<t>) - 101). (1) 

The covariant derivative 

D^ = d^ + ieA^ (2) 

gives a charge +e. We place only a few conditions on the Higgs potential V. 
There must be a minimum at zero for the unbroken- symmetry phase, and 
another minimum at some Vq for the broken phase. For a stable domain wall 
to exist we need V{vo) = ^(0); for convenience we take the common value 
to be zero. 

We expand the field 4> into real and imaginary parts as 

0=^(01+^2). (3) 

The Lagrangian becomes 
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Minimizing the action gives equations of motion for the three fields 0i, 02 
and A'': 



□01 = e<j)2d- A + 2eA^d^<P2 + e^A^<l)i-V'{J<l)l + <l)l) , (5) 

^1 + <PI 



□02 = -e0i9- A -2eyl^a^0i + 6^^202 -y'(V'02 + 02)^^4= (6) 



J2_ 
'I + 



d^F^" = e(0i9'^02 - 029'^0i)+eM'^(0? + 02). (7) 

These equations of motion have a stable domain wall solution with 0i = v{z) 
and 02 = = 0. We will find this solution and study small perturbations 
about it. We write 

01 = v{z) + h (8) 

and take h, 02 and to be perturbatively small. 

The term of order zero in equation (|^) is the condition for the field con- 
figuration v{z) to be stable: 

v" = V'{v). (9) 

We allow the prime symbol to have two meanings: since \^ is a function of v 
we take V to be dV/dv; in all other cases a prime indicates a derivative with 
respect to z. Integrating equation twice gives the domain wall solution 

^=r'"^i^ (10) 



in which the unbroken phase is taken to be at 2; = —00. We will also make 
use of the solutions v = Vq and v = which describe uniform broken and 
unbroken phases. 

The first-order terms in equations consist of an equation for the 

field h 

ah = -V"{v)h (11) 
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and the following coupled equations for 02 and A'^: 

V'iv) 

V 

d^F^" = e{vd^'(j)2-(p2d^v) + {evyA^'. (13) 



□02 = -ev d ■ A - 2eA>'d^v ^02 (12) 



If we take the divergence of equation (|I3|) and apply the condition (H) on the 
background field we get equation (|12D times a factor ev. Thus equation (|l^) 
is redundant as long as v is not zero. 

2.2 Particle content far from the wall 

To obtain the first-order equations appropriate to the phase of unbroken 
symmetry we set w = in equations (Pl)-(p!3|). Equation (|ll]) for h becomes 

{n+m^)h = 0, (14) 

where the mass m is defined by 

= V"{0). (15) 
Since v is zero, equation (|12[) is not redundant; it reduces to 

(□+m2)02 = O. (16) 



Equations ( p^ and ( [T6| ) combine to give a single Klein-Gordon equation for 
the complex field 0: 

(□ + m2)0 = O. (17) 



The final first-order equation ([T3|) simplifies to 

d^F^"" = 0. (18) 

These perturbative equations tell us that the unbroken phase contains a 
complex scalar field of mass m and a massless vector boson A'^. Our 
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linearized equations of motion describe only the interaction of particles with 
the background field f , not the interaction of particles with each other, so we 
must refer back to the original Lagrangian to find out that the scalar field 
carries charge +e. 

We define the following normalized solutions of equations (0) and (HI 
Their meaning and use are discussed below. 



0+ 






(p = ie^P-^ 


At 


A^' = ef'e-'P-'' 



(19) 



We choose our coordinate system so that the four- vector is {uj, p, 0, k). The 
energy u is always positive, but the momenta p and k may have either sign. 
The solutions (|T9|) are normalized so that their time-averaged 2;-component 
of energy flux is ^uk. In the case of the gauge boson solution At the nor- 
malization assumes that the polarization vector e'^ satisfies 



-1. 



(20) 



Although we use complex notation, the field A^^ is real-valued; the operation 
of taking the real part is implied. We subdivide At into two polarizations 
Ati and At2 with polarization vectors 



Eti 

^T2 



(0,fc,0,-p)/Vp2 + ^2 

(0,0,1,0). 



(21) 
(22) 



Using these standard solutions we can describe a field configuration with 
a single complex number. For example, the complex number Ati represents 
the solution 

A^ = ATis'^Tie-"'"'. (23) 



When a standard solution has negative energy we define the coefficient of 
the solution to be the complex conjugate of the given complex number. For 
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example, the complex number represents the solution 



(24) 



Next we consider the phase of broken symmetry. Setting v = Vq we find 
that the first-order equations (p!l|)-(]T3|) reduce to 

{n + m'^)h = 



(25) 
(26) 



where 



Defining the variable 



fh^ = V"{vo). 



ev 



(27) 



(28) 



which is invariant under perturbatively small gauge transformations allows 
us to simplify equation p6|) to 



where the mass M is given by 



(29) 



(30) 



According to equations (|25|) and (p9|) the broken phase contains a Higgs 
scalar h of mass m and a massive vector boson B'^ of mass M. 

The normalized solutions of the broken-phase equations are as follows: 



h 

B 



h 

B^' 



-ip-x 



(31) 



In addition to transverse modes Bti and Bt2 with the polarizations Eti 
and 5t2 given above, the massive vector boson has a longitudinal mode Bj^ 
with polarization 

'2 - M2 



UJ 



M2 



UJ 



-,p,0,/c). 



(32) 
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2.3 Separation of scalar internal modes 



We will reduce the equations of motion ([Tl|)-(p!3D in the domain wall 
background to independent scalar equations of the form 

[□ + U{z)] 5 = 0, (33) 

where S" is a scalar mode and U is the potential it sees. We take all fields 
to be eigenstates of energy and transverse momentum, that is, to contain a 
factor e~*'^*+*P^. Once we define the positive quantity E by 

E2 = cj2-p2^ (34) 

equation (^) becomes the Schrodinger equation 

[-dl + U{z)\ S = E^S (35) 

for a particle with nonrelativistic energy E'^, so our intuition for one-dimen- 
sional scattering can be applied to the potential U . Equation (|33|) shows that 
the asymptotic values f/(±oo) are to be identified with the mass squared of 
the scalar mode in the broken and unbroken phases. 

Equation (|ll]) for h is already in the required scalar form — we need only 
define the potential 

Un{z) = V'iviz)). (36) 



To extract scalar modes from equation (0) for A^, the first thing we do 



is write the equation in terms of the variable 5^ (28 



d^F^" = {evfB'^. (37) 

We can get two scalar modes by making the ansatz = 0. In this case, 
taking the divergence of equation (|37|) produces the relation d- B = 0, which 
allows us to simplify equation (|37|) to 



- nB^" = {evfB^". (38) 
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Factoring out a constant polarization by writing 



5^ = e'^r (39) 
gives a scalar equation for r, with potential 

Ur = {evf. (40) 
The constraints on B'^ allow two polarizations: 

Erl = {p,uj, 0,0) /E (41) 
Sr2 = (0,0,1,0). (42) 
We name the two corresponding scalar modes B^-i and 5^2- 



To obtain a third scalar mode from equation (37) we make the reason- 
able assumption that the polarization of the third mode is orthogonal to 
the polarizations of the two we already have. Referring back to equations 



(|T) and (1^) we find that orthogonality requires the first three components 
of the polarization of the third mode to be proportional to (cij,p, 0), so we 
can write 

B^" = d^'b for ^ 3. (43) 

Now the four components of equation (^7^ reduce to the following two equa- 
tions: 

[-dl + (evY] b = dsB^ (44) 
[-E^ + {evy]B'' = E%b. (45) 



Solving for b in terms of B^ yields^ 

1 



E^{ev) 



d,[{evrB']; (46) 



^We must eliminate b rather than B^, for the foUowing heuristic reason: Suppose 
we send a particle with the third polarization toward the wall with energy low enough 
that it will be totally reflected. When the particle reaches the classical turning point its 
momentum is {i-tJ,p, 0, 0). Now, for a massive vector boson with this momentum the third 
polarization is (0,0,0,1), so at this particular point in the scattering b will be forced to 
zero, not by the differential equation it obeys, but by the changing polarization vector. 

has no such singularity. 
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02 
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Figure 1: Connections between internal and asymptotic modes. 
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Bl 
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Bt2 



the remaining equation of motion is 

]ds + ievY + 2 



= 0. 



Defining 



B^ = {E/ev)\ 



removes the linear derivative term, leaving 



-E^ - dl + ievf + 2 



l\ 2 



V 
V 



A = 0, 



a scalar scattering equation with potential 



V J V 



(47) 
(48) 

(49) 
(50) 



This scalar mode, which describes the third polarization of the field B^, we 
refer to by the name Bx- 

2.4 Connection of internal and asymptotic modes 

A schematic diagram of how the various modes connect to one another is 
given in figure |l|. We will devote the rest of this section to making a precise 
statement of what this diagram means. 

The fields h and B^ which describe the scattering process are the same 
as the fields in the broken phase, so to find the connection between modes 
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Table 1: Summary of polarization vectors. The vector Ex is only meaningful 
in the broken phase. The last column contains normalization factors for the 
normalization = —1. 



at +00 we need only express the asymptotic polarizations of B'^ in terms of 
the internal ones. We summarize the various polarization vectors in table |l[ 
The only new information in this table is the polarization vector Ex. To 
calculate this vector, we consider a scalar function A which is asymptotically 
a plane wave: 

A = e-'P-". (51) 
Working backward through equations (^Sj) , (^6|), and (^) we find 

53 = (l)e--- (52) 

b = -J—(ik + -)e-'P-'' (53) 
Eev)\ vj ^ ^ 



P^ / y 



B^" = ^--{k-i-]e-'P-\ (54) 
E{ev)\ vJ ^ ^ 

Using the fact that v' /v — > in the broken phase, we read off the unnormal- 
ized polarization vector from equations ( p^ and (|5^). 

Expanding the broken-phase polarization vectors in terms of the internal 
ones we obtain the following matrix equation connecting the internal modes 
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to the broken-phase modes: 

( h \ /I 

Bl 



Bti 

V Bt2 / 



pM 



V 



pM Lok 
Ei/p^+fca E^p'^+k'^ 



\ I ^ \ 

Bx 



1 / 



BtI 
V Br2 J 



(55) 



where zeroes in the matrix have been replaced by blanks for legibility. At 
normal incidence the matrix reduces to the identity. 

The connection between the internal modes and the asymptotic modes in 
the unbroken phase is more involved. We begin with the simplest case, the 
connection between Br2 and At2- We consider a solution of the scattering 
equation for the 5^2 mode which is asymptotically a plane wave, that is, for 
which there is a complex number Br2 such that 



r - Br2e 



-ip-x 







as 



z — > — oo. 



(56) 



According to the definition (|39| ) of r, this solution generates a field configu- 
ration 

B^[Br2]=e>',^T. (57) 

On the other hand, from the definition (p8|) of B^ we find that the field 
configuration generated by the standard solution for the At2 mode is 



T2 



^T2 ^T2^ 



(58) 



If we take At2 equal to i?^2; the two field configurations are asymptotically 



the same, that is. 



B''[B^^] - B^'[At2] ^ Q as 



-oo. 



(59) 



We use this condition — that the difference between field configurations go to 
zero — as our criterion for connection. 
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We handle the Bri mode similarly. An internal solution with asymptotic 
form 

T = Brie^'P-"" (60) 
generates an asymptotic field configuration 

B'^[Br^] = B,,e^,,e-'P■^. (61) 

Applying the relation 

e^^i = (sgnA;)4,+ (-^y (62) 

separates B^[Bri] into two terms. The first term matches the field 
generated by an Ati solution with amplitude Ati = {sgnk)Bri. The second 
term describes a field 

= Bri {£^P^^~^""' (63) 

polarized in the direction of p^. Although the field ( |63D is a valid solu- 
tion of the unbroken-phase equation of motion (|THp for A^, it has a field 
strength {F^^) of zero and so carries no energy or momentum. We discard 
this pure gauge solution and find that B^-i connects directly to Ati- 

The use of the coefficients B^-i and Bt-2 in the preceding paragraphs has 
obeyed the following convention, which we will continue to observe: whenever 
we use a complex number to describe the asymptotic form of one of the 
internal modes, the number is the coefficient of a unit-amplitude plane wave. 
We will apply the same convention to h and 02 even though they are not 
internal modes. 

Next we discuss the connection between the internal mode Bx and the as- 
ymptotic modes in the unbroken phase. A solution of the scattering equation 
which is asymptotically a plane wave, that is, for which 

A - Eac-^*"^ ^0 as -oo, (64) 
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generates an asymptotic field configuration 



B'[Bx]=B,i-)e-'P-''. (65) 



.evj 

Since w ^ in the unbroken phase, this field configuration is divergent. 
However, we can obtain an equal divergence from the asymptotic field 02- We 
take 02 to be a plane wave with amplitude given by the complex number 02 
and substitute into the definition (|28|) of -B^ to find 

B-'[<P,] = <pJ-){--tk)e-^P-^. (66) 
\evj V 

As f ^ 0, the ratio v'/v reduces to m, so we can match the two divergences 
by taking 

• ^5 .B, = '^B,. (67) 



m — ik E 

Although equation ( |67|) will turn out to be the correct relation between 
the modes Bx and 02, the above argument using matching divergences does 
not prove that the two modes connect, because the difference of the matching 
divergences contains a finite part which does not necessarily go to zero as 2; — *■ 
—00. To calculate the finite part of the field B^^lBx] correctly, we keep the 
first-order term in the asymptotic expansion of A near v = 0, writing 

X = Bxil + cv)e-'P-'', (68) 

where the constant c is determined by the scattering equation (|49|) for the 
Bx mode. A calculation which we omit shows that when the relation between 
the modes Bx and 02 is that given in equation (|67D, the field configurations 
Bi' lBx] and 5'^[02] satisfy 

B^'[Bx] - (fi^[02] + Cp^e-'P-"") ^0 as z ^ 0, (69) 

where the coefficient C of the pure gauge solution is an appropriately chosen 
constant. In other words, the field configuration B^^lBx] is asymptotically 
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equal to the field -B'^[02] plus a pure gauge term. We discard the pure gauge 
solution and find that Bx connects to 02 as in equation (p7|). 

All that remains is to combine the 02 produced by Bx with the internal 
mode h to get the unbroken-phase modes 0^ and 0~. Since h and 02 are 
real-valued fields, the amplitudes h and 02 actually represent the functions 



and 

which we combine into 

0[/l,02] = 



h 


-ip-x _j_ 




— e 


— e 


2 


2 


02 - 






— e 


-ip-x _j_ 


% 


2 


2 



ip-x 



h + 
~V2 



,1 . 



+ 



h* + z<; 



:le 



(70) 
(71) 

(72) 
(73) 



Using the normalized solutions (|l^) we read off the amplitudes 

. h + i(j)2 , A- ^~ ^<^2 

We collect all our results about the connection between the internal and 
unbroken-phase modes in the following equation: 

/ 0+ \ / 1/V2 \ f h \ 

1 / /o k—im 

sgn k 

1 y 



V 



5a 

V 5-2 y 



(74) 



Ati 

V ^T2 / 

The matrix equations (|55D and (|74D are the desired precise statement of what 
the connection diagram in figure |I] means. 



3 Scattering from the domain wall 

3.1 Specification of the potential 

All our results up to this point are independent of the details of the Higgs 
potential V{v). To obtain anything more than generalities about scattering 
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from the domain wall, however, it is necessary to give up this independence 
and specify a form for the potential. We choose the standard quartic potential 

1 



V{v) = -\v\v-vof. 



(75) 



In section ^.7| we will discuss how our results depend on this choice. 

With this choice of potential the mass m ([T^) of the charged scalars is 
equal to the mass m (|27|) of the Higgs, 



2 ~ 2 

m = m 



0; 



(76) 



and the integral (|T0D which determines the background field can be evaluated 
explicitly, with result 



v{z) 



1 + tanh 



mz 



(77) 



The scattering potentials Uh, and U\ can be cast into the standard form 



Uo{l - s) + UiS + U2S{1 



where the dimensionless coordinate s is defined by 



mz 



-I 1 + tanh — I = — . 

^^0 



(78) 



(79) 



The coefficients Uq, Ui and U2 for the various potentials are listed in table |[ 



3.2 Solution using hyper geometric functions 

In this section we solve the scattering equation 

[□ + u{z)] s=o m 

analytically for any potential U of the standard form ([78|), that is, any po- 
tential quadratic in the variable s. Although our results are more general. 



15 





General V 


Specific V 


Uo f/i t/2 


h 

Br 

Bx 


Uh = V"{v) 
Ur = {evy 

Ux = {evr + 2{^-y-^i 


m2[l-6s(l-s)] 
m2(l - s) + M^s^ 


1 1 -6 

-M^ 

1 -M^ 



Table 2: Summary of internal mode scattering potentials. The mass ratio Ai^ 
is defined to be M'^ /m?. 



our method of solution is similar to that used by Ayala et al to study the 
scattering of fermions. We factor S into 



S = s'^'^l- sY^ S{s)e 



—iu)t+ipx 



(80) 



where kq and Ki are constants to be determined, and find that the scattering 
equation becomes 



2kq s{l - sf - 2ki s^{l -s) + s{l - - 2s; 

+ ( Kq^ (1 - sf + Ki' - (kq + Ki+ 2KoKi) s(1 - s) 

+ — -Uo{l-s)-Uis-U2s{l- s) 



d 
ds 



Setting 



l^i — Ui - 



S = 0. (81) 



(82) 



makes equation (|8TD divisible by s(l — s); the quotient is the hypergeometric 
equation^ 



s(l -s)— + (^7- (a + /3 + l)s^^ - a/3 



S = 



(83) 



^AU the relevant facts about the hypergeometric equation appear in reference [^. 
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with parameters 

a = kq + Ki + ^ + 5 



P = Ko + /ti + |- (5 (84) 
7 = 1 + 2ko. 



The constant 6 is defined by 



6 = Jl/4 - U2. (85) 



One solution of the hypergeometric equation (|83| ) is the hypergeomet- 
ric function F(a,/3, 7;s), and the corresponding solution of the scattering 
equation (^) is 

S = (1 - s)^' F{a, P, 7; s) e-^"*+^f^ (86) 

We study the asymptotic behavior of this solution, beginning with the case 
z — > — 00, that is, s ^ 0. As its argument s goes to zero the function F 
reduces to unity, as does the factor (1 — 5)*^^ From equation ( [7^ ) we find 
that as 2; — > —00 the coordinate s simplifies to e*"^, so the solution (|86D is 
asymptotically the plane wave 

^—iojt+ipx+mKoz (87) 

For the solution to describe scattering, this exponential must represent the 
transmitted part of a wave incident from z = +00, that is, kq must be 
negative imaginary. To avoid confusion, we take Ki to be negative imaginary 
as well. 

To calculate the asymptotic behavior of the solution ( p6D in the opposite 
limit z +00 we use the identity 

F{a, P, 7; s) = r(7)r(7-«-^) „ + ^ + 1 _ ^. 1 _ (88) 

r(7 - a)r(7 - /?) 

+ ^^^^r{a)m~ ~ ^)'"""'-^(^ -a,^-P,^ + l-a-P;l-s). 



17 



Taking into account the parameter values (p^ ) and replacing 1 — s by its 



equivalent e ™^ we find S has the following asymptotic form at +00: 



T(1 + 2ko)T(-2ki) 



T{ko - Ki + ^ + 6) r(Ko - Kl + I - 5) 



r(i + 2Ko)r(2Ki) . ^ 

Since ki is negative imaginary, the second term represents the incident wave. 
Taking the appropriate ratios yields the refiection and transmission coeffi- 
cients 

T{ko + K1 + I + 5) T{ko + + r(-2Ki) 



r(/€o-/ti + | + 5)r(/€o-/ti + |-5) r( 2/ti) 



(90) 



^ ^ r(/to + /«i + | + 5)r(/to + fi:i + |-5) 
r(i + 2/to) r(2/ti) 



In the case of total refiection we make the transmitted wave (|87|) into a 
decaying exponential by taking to be positive real. Although the particle 
is totally refiected, the transmission coefficient t does not become zero, and 
indeed it need not, being merely the coefficient of a decaying exponential. 
The refiection coefficient r does reach unit amplitude, however. 

The solution of the hypergeometric equation (^) which represents a par- 
ticle incident from — cxd is 

S = F{a,p,a + p + l-rA- s). (92) 

Proceeding as before we find refiection and transmission coefficients which 
differ from the previous results (pO|) and (|9T|) only by exchange of Kq and Ki. 

3.3 Bound states 

To search for bound states we take both hq and ni to be positive real and 
require that the growing exponential in equation ( pUD have coefficient zero: 

r(i + 2/to)r(2Ki) 



0. (93) 
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The gamma function has no zeroes, but it does have a pole at each nonpositive 
integer. The parameter a = + + ^ + 6 is always either strictly positive or 
complex, so equation ( PBD reduces to the condition that the parameter (3 = 
kq + Ki + ^ — 6 he a nonpositive integer, that is, that 

- P = 6 - I- {ko + Ki) = n (94) 

for some nonnegative integer n. We attempt to satify this condition by 
varying the energy E'^. Our variations must keep kq and Hi real- valued and 
strictly positive, which, according to the definition (|8^), means that 

— < min (f/o, f/i) (95) 

and therefore that 



Ko + Ki> ^\Uo-Ui\. (96) 

As a result, —f3 may take any value from — oo up to but not including a 
maximum given by 

N = 6-l- ,J\Uo-U,\, (97) 



and equation (Q) can be satisfied for any nonnegative integer n < N. Solving 
equation (0) for E"^ yields the energy of the corresponding bound state 

For bound states the solution (|86|) of the hypergeometric equation simplifies, 
because when — /3 is equal to a nonnegative integer n the power series for the 
hypergeometric function 

7 7(7 + 1) 2! 

collapses to a polynomial of degree n. The lowest-energy bound state has 
the particularly simple wavefunction 

^ = s"«(l - s)"^ e-^-^^+^P^ (100) 
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Using the parameters Ui given in table we can compute the number of 
bound states for each of the internal modes. For the mode we find the 
bound state parameter is never positive, so there are no bound states, but 
then we expect none in a potential without a minimum. 

The h mode has N = 2 and so has two bound states. We apply equa- 
tions (|98|)- (|100|) to find that the first bound state has energy = and 
wavef unction 

h = s{\- s) e-^-*+^P- = - sech' ^ ^-Mpx ^^q^^ 
and that the second has energy E'^ = (3/4)m^ and wavefunction 



/ 1 TT) 7 7117 

h = ^s{l - s) (1 - 2s) 6-*"*+**^" = -- sech tanh e-'^'+'P\ (102) 

To make a physical interpretation of these bound states we remember that 
h represents a small perturbation added to the background field v. The fac- 
tor sech^(m2;/2) which appears in the first bound state ( |101| ) is proportional 
to v'{z)^ and adding a small amount of v'{z) to v has the effect of translating 
the wall in the 2;-direction, so the wavefunction ( |101| ) describes a plane-wave 
oscillation of the local position of the wall. The second bound state ( |102| ) 
describes a similar oscillation of the local thickness of the wall. These two 
excitations of the wall can be thought of as quasiparticles confined to the 
wall surface with masses u"^ — = E"^ of zero and (3/4)m^ respectively. No 
matter what Higgs potential V we choose, the h mode always has a bound 
state for oscillations of the wall position — we need only apply the equation 
of motion (P) for v to find that the wavefunction h = y'(^z)e~'^'^'^~^^^^ satis- 
fies equation (|lTD whenever E^ = 0. The underlying physical reason for the 
existence of this bound state is the translation invariance of the Lagrangian. 

For mass ratios Ai^ = M"^ /m? in the range 3/4<A^^<2 the B\ mode 
has a single bound state. Neither the energy ( pHD nor the wavefunction ( |100| ) 
of this state simplify further. The bound state parameter N and energy E"^ 
are plotted as functions of in figure H. 
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— ^ ^ \ ^ 1 1 ^ ^ ^ l-M"^ 

1 2 3 

Figure 2: Bound state parameter N and bound state energy E"^ /w? as func- 
tions of Ai^ for the B\ mode. The dotted hne in the second plot marks the 
start of the continuum of unbound states. 
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3.4 Reflection and transmission probabilities 

For the case in which total reflection does not occur we calculate the reflec- 
tion and transmission probabilities R and T for a particle incident from +00. 
Starting from the coefficients (|90|) and ( pi]) and applying the identity 



we find 



Viz)V{l-z) = -!^ (103) 
smTT^; 



2 sin7r(Ko - Ki + i + (5) sin7r(/to - Ki + I - 5) , 

R=\r\^ = ^ ^ (104) 

sin 7r(Ko + + 2 + ^) sin 7i{ko + Ki + ^ ~ ^) 

ko . 12 — sin 27rKo sin 27rfi:i 

~ ki ~ sin7r(Ko + Ki + i + 5) sin7r(Ko + + I - 5) ' 

where the asymptotic momenta ki are given by 



^2 _ jYi2lf. — ifYii^._ (106) 



The probabilities ( p.04| ) and ( |105D are invariant under exchange of kq and ki 
and so apply to particles incident from —00 as well. 

Using trigonometric identities we convert the transmission probability 
from equation (|105|) into the more practical form 

2 sinh ^ sinh ^ 
cosh ^ cosh ^ + sinh ^ sinh ^ + cos 27r5' ^ ^ 

mm mm 

Except in a few cases which we discuss momentarily, the transmission prob- 
ability is a strictly increasing function of the energy E"^ with a square-root 
singularity at the threshold of total refiection. A typical curve is shown in 
figure 1^. When Uq = Ui the square-root singularity is replaced by a linear 
approach to zero, unless also cos27r5 = —1, in which case the transmission 
probability T = 1 is no longer even strictly increasing. According to figure ^ 
the energy E"^ need only be a small fraction of the mass above threshold for 
the transmission probability to be significant. Other choices of the scattering 
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Figure 3: Transmission probability for the B^^ mode when Ai"^ = 1. 

potential do not make the onset of transmission significantly slower — in all 
cases the transmission probability reaches 1/2 before E"^ gets 0.02 above 
the threshold. On the other hand, the onset becomes much faster as the 
potential approaches any of the potentials which give T = 1. 

The transmission curves for the Br and B\ modes have the generic form 
described in the previous paragraph. If the gauge boson mass M were zero, 
the Br mode would have T = 1, so if the gauge boson mass M is small 
compared to m the curve for the Br mode will have a sharp onset. The 
onset in the Bx mode is never particularly sharp. For the h mode, Uq = Ui 
and 5 = 5/2 make the transmission probability unity. This result does not 
generalize to an arbitrary Higgs potential V. For most choices of V, the 
mass m of the charged scalars is different from the mass fh of the Higgs, so 
that reflection — in fact, total reflection — is possible. 
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3.5 Scattering of asymptotic modes 



In section p.4| we derived the connection equations 



/ h \ 

Bl 

Bti 

V Bt2 



( 1 
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pM 
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pM uik 



E^p^+k'^ E^p^+k'^ 



and 
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V j 
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k—im 
' Eyjl 
k—im 
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( h \ 
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V 5r2 J 



sgn k 
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h \ 
Bx 

Brl 
V Br2 J 



m 



Using these two equations and the reflection and transmission coefficients 
( pOD and ( pT] ) we can calculate scattering probabilities for a particle incident 
from the broken or unbroken phase. We use an inverted connection equation 
to convert a unit amplitude for the incident asymptotic mode into amplitudes 
for internal modes, multiply by the appropriate reflection and transmission 
coefficients to find the scattered amplitudes at both positive and negative 
infinity, and use the connection equations to convert back to amplitudes for 
asymptotic modes. We take the norm of these amplitudes to get probabilities, 
multiplying in the case of transmission by the appropriate momentum ratio. 

The meaning of the symbol k depends on the context in which it appears. 
To avoid confusion, we define for each of the different particle masses a 
corresponding positive momentum, for example. 



;io8) 



and write each k in terms of these momenta. In the connection equation (^) 
for the broken phase, k is the momentum of a massive vector boson, that 
is, \k\ = ku- Because equation (|55| ) expresses asymptotic amplitudes in the 
broken phase in terms of internal ones, we will apply it to particles moving 
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away from the wall, that is, particles with positive 2;-momentum. Setting 
k = +kM we obtain 



/ h \ 

Bl 

Bti 

V Bt2 J 



( 1 



V 



ujkM 



pM 



pM 



\ f h \ 

Bx 



1 / 



BtI 
V 5r2 J 



(109) 



To obtain the connection equation for particles moving toward the wall, we 



invert equation (|55D and take k = —ku, with result 
/ 1 



/ h \ 

Bx 

BtI 
V Bt2 ) 

Substituting k 



u)kM 



pM_ 



pM u)kM 



\ f h \ 

Bl 



V 



1 / 



Bti 

V Bt2 



(110) 



-km into equation (|74D gives the connection equation for 
particles moving away from the wall into the unbroken phase, 

/ (J)+ \ ( 1/v^ \ ( h \ 

T I ^f^ km+im 

-1 

1 J 



V 



Ati 

V ^T2 J 

while inverting and substituting k 
incident from the unbroken phase: 



Bx 

BtI 
V Bt2 J 



(111) 



+km gives the equation for particles 
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1/V2 1/V2 

km+im km+im 
EV2 EV2 



\ 



At2 J 



(112) 



We attach subscripts to the symbols r, t, and T to indicate to which 
internal mode they apply. As before, the symbols r and t represent reflection 
and transmission coefficients for a particle incident from +oo, so for example 
we have 



h 

kjM 



ty 



(113) 
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Using the method described at the beginning of this section, we calcu- 
late the transmission probabilities for a Bl particle incident from the broken 
phase. We follow the initial amplitude vector through the processes of con- 
version into internal modes, transmission, and conversion into asymptotic 
modes at — oo, finding 



/ \ 
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V / 
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(111) 
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To obtain transmission probabilities we take the norm of each element of 
the final vector and multiply each by the appropriate momentum ratio. The 
resulting transmission probability vector is 



_1 

2 kM 
2 kM 
kM 
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\tr\ 
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;ii5) 



where the probabilities P and Q are 

P = 

Q=l-P = 



Mj 



(116) 
(117) 



MJ 



The normal connection probability P ranges from unity at normal incidence 
to zero at grazing incidence; its behavior at intermediate angles is shown in 
figure ||. 
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Figure 4: Contour plot of the normal connection probability P as a function 
of p/M and kM/M, with contour interval 0.1. The probability P is unity at 
normal incidence and zero at grazing incidence. 

The results ( [115|) of the previous calculation suggest a simpler method 
which we describe for the case of transmission in the 0"'" mode. As indicated 
in figure ^ there is only one path from Bl to 0"^. Since there is no interfer- 
ence between different paths we can calculate with probabilities instead of 
probability amplitudes and can read off the transmission probability PT\/2 
directly. Using this method we obtain most of the scattering probabilities 
given in table ^. For generality we do not make use of the fact that Th = l. 
The only entries of table ^ that require further calculation are the ones that 
describe reflection of the modes and Bti and reflection of the charged 
scalars </>"'" and . 
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Figure 5: The unique path from Bl to 0"*". 
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Table 3: Scattering probabilities for asymptotic modes. A blank entry indi- 
cates a probability of zero. 
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3.6 Interference between internal modes 



We consider the reflection of an incident particle. Tfie amplitude 
vectors for the process are 



/ \ 
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V J 
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(110) 
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ujkM 



E^P^+kli 
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(118) 
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We take the norm of the elements of the final amplitude vector to obtain two 
reflection probabilities, the probability -R^ of reflection of 5^ as itself and 
the probability X of crossing, that is, of reflection of into -Bti- When 
we calculate the corresponding reflection probabilities for a Bti particle, we 
find that the crossing probability X also describes reflection of B^i into B^. 
The reflection probabilities are 



Rl 
Rti 
X 



P^Rx+ Q^Rr-2PQ^RxRrX 
Q^Rx+ P2 Rr-2PQjRxRrX 



PQ Rx + PQRr + 2PQjRxRr X 



where the crossing enhancement factor X is defined by 

Re 



X 



X'r 



\rx\\rT 



(119) 
(120) 
(121) 

(122) 



The interference term which appears in equations (|119D -( p!2T|) is pro- 
portional to the factor X, the cosine of the angle between the amplitudes 
Tx and Tt-. When X is positive, the interference term increases the proba- 
bility of crossing. For mass ratios A^^ > 1 the behavior of X as a function 
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Figure 6: Crossing enhancement factor X as a function of E'^/m'^ for a mass 
ratio M'^ of 2. 



of E"^ /m^ is not complicated. The curve shown in figure |^ is typical — X drops 
from the value 1 at the threshold E"^ /m? = M.^ to a positive minimum, then 
returns gradually to 1. 

We plot a set of representative curves for mass ratios < 1 in figure |^ 
In the region E"^ /m? > 1 the curves all exhibit the same positive minimum 
and gradual return to 1, but in the region E"^ /m? < 1 where the Bx mode is 
totally reflected, the curves are qualitatively different from one another. The 
mass ratio A^^ at which the cusp at E'^ /m? = 1 reaches unit height is only 
approximately 0.4225, but the ratio at which the X = —1 minimum vanishes 
is exactly 3/4, the same ratio at which the bound state in the Bx mode 
appears. 

Analogous calculations for the reflection of the charged scalars and 0~ 
yield the reflection probabilities 



2\ RhR 



\[Rh 

\ {Rh + Rx + 2^ RhR\ X(j, 



where 



X^ = Re 



E 



Fh rx 



(123) 
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Figure 7: Crossing enhancement factor X as a function of E'^/m^ for gauge 
boson mass ratios M"^ of 0, 0.4225, 0.7, 3/4, 0.8, and 1. 



For the quartic Higgs potential (|75D these resuhs for charged scalars simph- 
fy — the amphtude Vh is identically zero, X</, is undefined, and the reflection 
probabilities and both reduce to -Ra/4. 

3.7 Discussion of results 

For particles with energy sufficiently far above all relevant mass thresh- 
olds the internal mode transmission probabilities are essentially unity and 
the process of scattering from the domain wall can be regarded as a pro- 
cess of conversion of particles between phases. In this limit the scattering 
probabilities in table ^ reduce to the conversion probabilities in table ^. The 
breakdown of this approximation occurs when either the momentum of the 
incident particle approaches zero or one of the internal modes to which the 
particle connects approaches total reflection. 

We consider the extent to which our results based on the quartic poten- 
tial of equation ( |75D apply when the Higgs potential is more general. With 



the exception of the plots of X, the results of sections IXB and PTB] hold for 
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Table 4: Conversion probabilities for high-energy particles. 

an arbitrary Higgs potential — all dependence on the potential is contained 
in the reflection and transmission amplitudes for the internal modes. From 
the discussion of bound states in section |3]^ we draw the following generally 
applicable conclusions: The potential Ur is always monotone increasing, so 
the Br modes never have bound states. The Bx mode need not have a bound 
state, but is capable of supporting at least one. The h mode is guaranteed 
to have at least one bound state. We cannot make any more specific state- 
ment about the bound states of the h mode, because for any positive integer 
there exist Higgs potentials V{v) for which the h mode has that number of 
bound states. Suitable Higgs potentials can be constructed by joining three 
quadratic pieces; the resulting scattering potentials Uh = V"{v) are square 
wells of various sizes. 

4 Extension to the standard model 

4.1 Reduction to two subproblems 

The fields appearing in the standard model are listed in table ||, along 
with their transformation properties under the various symmetry groups. We 
take the symbols B^ and, later, A^^ to have completely different meanings 
than in the previous sections. 
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Table 5: Fields appearing in the standard model. A check indicates that 
the field transforms according to the fundamental representation of the cor- 
responding group. 



The standard model Lagrangian is 
L = 

4 

- 4Me0e^ - q^Mu(t)cUn - q^Md(j)d^ + H. c. ). 
For equation ( p.26| ) we define the norm |0| by 



(126) 



101 = V^ 



and the conjugate field 0c by 



(127) 



:i28) 



The constants Mg, M^, and are matrices in family space. 

To study scattering from a domain wall at the electroweak phase transi- 
tion, we will write out the components of 0, 

,^ ^1^03+^4^ 
V2V01+#2. 



(129) 
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derive equations of motion, define 

take all fields except v to be perturbatively small, and obtain an equation 
for the static domain wall solution v along with first-order equations of mo- 
tion which describe perturbations about it. First, however, we simplify the 
calculation by removing from the Lagrangian terms which do not contribute 
to the end result, such as terms higher than quadratic in fields other than 0. 

Specifically, we remove the bilinear terms from the gauge field strengths 
and the gauge couplings from the fermion covariant derivatives. The fermion 
mass terms are quadratic in the perturbatively small fermion fields and so 
do not contribute to the equations of motion for cf). Consequently, we can 
expand according to equations ( |129| ) and even before deriving the equa- 
tions of motion. Dropping higher than quadratic terms yields the simplified 
mass terms 

-y|-efl -Mi-^Mfl + H. c.j. (130) 

Redefining the fermion fields by applying independent unitary rotations in 
family space to the left and right components of the fields e, u, and d allows 
us to take the mass matrices Me, M^, and to be diagonal with positive 
real eigenvalues. 

As a result of the above transformations the Lagrangian (|126D falls into 
unrelated pieces. For each massive fermion species there is a copy of the 
Dirac Lagrangian 

^ ^ + H. c.) - {fxv)i;ip (131) 



.2 

with a broken-phase mass hvq equal to times the eigenvalue of the 

appropriate mass matrix. We consider the massive fermions further in sec- 



tion |4.3| . The pieces of the Lagrangian which describe neutrinos 

^uju, + H. c.) (132) 
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and gluons 

1 



{d^Gl-d^G''^){d^'G'"' -d^G""") (133) 



4 

are Lagrangians of free particles. To first order the neutrinos and gluons do 
not interact with the domain wall, and we do not consider them further. The 
remainder of the original Lagrangian, 

+ {D,<py{D^<j>)-V{V2\<j>\), (134) 
we analyze in the following section. 

4.2 Standard model bosons 



Having reduced the standard model Lagrangian ( 126 ) to the boson La- 
grangian ( |134|) , we apply the method described in the previous section and 
derive equations of motion. After some calculation, we find that the condition 
for the field configuration v{z) to be stable 

v" = V'{v) d) 

and the first-order equation of motion for the field h 

ah = -v"{v)h (0) 

are the same as in the Abelian Higgs model. The equations of motion for the 
fields 03 and W^", 

□ 03 = -^vd -W^ - gW^^'d^.v -^^^-^(Ps (135) 

where 

W^f^'^ = Qf^W''' - d^W'^", (137) 
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are exact analogues of equations (|T2D and ( ]T3| ) for 02 and A^, as are the 
equations of motion for 04 and W^^: 

V'{v) ^ 



□04 = -Igvd-W^ - gW^^'d.v - 

^ V 

d^W^^"" = lg{vd''(j)4-(l)4d''v) + {\gvyw^^. 



To cast the last three equations of motion 



□ 02 = -lvd-ig'B-gW'')-ig'B'^-gW^^')d^v 
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(139) 

02 (140) 
(141) 
(142) 



into a form analogous to equations (12) and (13), we define the linear com- 
binations 



g'B^'-gW^^' , ^„ gB^' + g'W^^' 



- - rr and 

^fg^TT" 

and obtain equations of motion for 02 and 



□ 02 



\Jg^ + g'^vd-Z- Jg^ + g'^ Z^^d.v - 



V\v) 



(143) 



(144) 



= + ^'^ (^^'^'^s - 02^^^^) + dV^?' + 9''' vf Z^ (145) 

along with an equation of motion 

d^F^"" = (146) 

which describes a noninteracting massless gauge boson — the photon. 

We summarize the correspondence between the standard model and the 
Abelian Higgs model in the following table. 



A^' 02 e 


B^' M 


W'^^ 04 y 

02 Wg-' + g'-' 
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(147) 
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The right side of the table shows the names we assign to the analogues of B'^ 
and M. From the defining equations (^) and (|30| ) for and M we find, 
for example, 



= W^f' + d^i^] (148 



2/i _ 

Ml = {Id'^of- (149) 



By applying the above correspondence to our results for the Abelian Higgs 
model we can find out nearly everything we want to know about standard 
model bosons. To be specific, we obtain the definitions and scattering po- 
tentials of the scalar internal modes and most of the equations that connect 
the internal modes to the asymptotic ones. What the correspondence fails to 
provide is the connection between the real Higgs modes 02? 03 and 04 and 
the complex modes appropriate to the unbroken phase. The modes h and 02 
combine just as in the Abelian Higgs model — copying equation ([73| ) we define 
the amplitudes of the particle and antiparticle solutions 

0° = ^^ and 00 = (150) 

and obtain the connection diagram given in figure ||. The associated connec- 
tion equations, scattering probabilities, and conversion probabilities may be 
found by applying the correspondence to equations (|55|) and ([7^) and tables 
I and |. 

Since the Higgs fields 03 and 04 combine with each other rather than with 
the field h in the charged particle and antiparticle amplitudes 

- and 0- = ^i^, (151) 



V2 72 

the connection equations involving these modes will not be exact analogues 
of the connection equations of the Abelian Higgs model. Taking the relevant 
parts of equations (|55| ) and (|73|) and applying the correspondence, we obtain 
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Figure 8: Connections between neutral modes in the standard model. 
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along with a similar pair for the fields W"^ and 04. We define charged gauge 
boson amplitudes for each of the internal and asymptotic gauge boson modes, 
for example, 



Tl 



Tl 



and 



Since the definitions ( |151| ) and ( p. 54] ) describe the same linear transformation 
of amplitudes, the connection matrices for the charged bosons are identical 
to those in equations (|152|) and (|153|) . We present one of the associated 
connection diagrams in figure § the absence of connections between positive 
and negative bosons is of course guaranteed by the conservation of electric 
charge.0 Applying the methods of sections |3.5| and produces the scattering 



''The correct analogue of the charge non-conservation in scattering in the Abehan Higgs 
model is the non-conservation of isospin and hypercharge, for example, in the reflection 
process 
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Figure 9: Connections between charged modes in the standard model. 
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Table 6: Scattering probabilities for charged modes in the standard model. 

probabilities given in table ^ and the conversion probabilities, that is, the 
transmission probabilities for high-energy particles, given in table |^. 

4.3 Scattering of fermions 



Variation of the Dirac Lagrangian (|131|) gives the Dirac equation 

{i$ - fxv)ip = 0. (155) 
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Table 7: Conversion probabilities for charged modes in the standard model. 
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To reduce this equation to scattering equations for scalar internal modes we 
use the ansatz of Ayala et al 

ip = {i$ + iiv)^p (156) 

to obtain 

[□ + {^ivf - il'^{^lv')\ = 0. (157) 

We take ip to be proportional to an eigenvector u± of the matrix — with 
eigenvalue ±1, 

i' = u±f±, (158) 

and find that the corresponding scalar function f± satisfies the scattering 
equation 

[□ + f//±(^)]/± = (159) 

with potential 

Uf± = {^ivf ±liv'. (160) 

The existence of two different potentials for fermion scattering is an ar- 
tifact of the process of extracting scalar modes, because any solution of 
equation (|155| ) can be expressed in terms of either /+ or For example, 
the solution ip generated by an eigenvector and function /_ can also be 
obtained from the eigenvector u+ and function /+ given by 




/+ = ^[^iv + d,y^. (162) 

We will express our results in terms of u+ and /+, because of the following 
advantage of the potential being the sum of two positive terms which 
go to zero at —00, it has no absolute minimum and so clearly has no bound 
states. 
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Although we have reduced the Dirac equation ( |155| ) to a scalar equa- 



tion ( p.59| ), we have not yet defined internal and asymptotic modes for the 
fermion field. For the moment we consider only positive-energy solutions of 
equation (|155|) . In the Pauli-Dirac representation of the matrices 7^, the 
normalized asymptotic solutions for the field ip can be written in terms of a 
two-component spinor (p as 

*i=^^f"+';")^.— . (163) 

Our normalization — setting the z-component of energy fiux equal to ^uk — 
is compatible with massless particles, so these asymptotic solutions can be 
used in both the broken and unbroken phases. The two-component spinor ip 
describes the spin of the fermion in its rest frame; our normalization assumes 
it satisfies 

(/?V = 1- (164) 

To study the relation between internal and asymptotic solutions we con- 
sider a solution which is asymptotically a unit-amplitude plane wave: 

/+ = e-*^-". (165) 

We write the eigenvector m+ in the form 

«+ = (1-^7')(q), (166) 

where X is a constant two-component spinor. Substituting into the definitions 
( |156D and ( |158D , we find that the resulting asymptotic field ip is equal to 



for the following unnormalized spinor (f: 

V. = ^i=f ^ + ^^-^^ , ]x. (167) 

As before, we choose the coordinate system so that the transverse momen- 
tum p lies along the x-axis. Because the values of the quantities k and v 
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appearing in equation ( |167| ) depend on whether the particle is in the broken 
or unbroken phase and on what the sign of the 2;-momentum is, the spin 
direction determined by (p can change as a result of the scattering process. 
However, the particular choices 



X 



generate spinors 



±i 



{u + fiv ^ p — ik) X 



:i68) 



(169) 



which point in the directions y and —y independent of the values of k and v. 
Accordingly, we define the internal modes /+| and f+i to have spins aligned 
with the y-axis, that is, to have the normalized polarization spinors 



2J2{u;-p) 



. and 



X 



2J2{lu + p) V2 



4. \X im 



Although equation ( |169| ) contains the nonconstant quantities k and v, the 
normalized spinors X^ and Xi are constant, as required, as a result of the 
identity 

\uj + fiv ^ p — ik\'^ = 2{lj + fiv){u =F p). 



Defining the asymptotic modes 



^1 =ip 



71' 



and 



71' 



(171) 



(172) 



we obtain the connection equation 




/ 



V 



Lj+fiv—p—ik 



Lu+fiv+p—ik 




(173) 



which is valid in both the broken and unbroken phases when the appropriate 
values of k and v are used. Because the connection matrix is diagonal, we 
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conclude that our chosen asymptotic modes do not interconnect and that 
the scattering probabihties for each of the asymptotic modes are equal to the 
scattering probabilities -R/+ and Ty+ for the internal mode /+. The appear- 
ance of nontrivial complex phases in the connection matrix of equation ( |173| ) 
indicates that in general the spin of the fermion rotates about the y-axis 
during scattering. 

To obtain the negative-energy solutions of the Dirac equation ( |155| ) we 
make use of the charge conjugation symmetry^ 



(174) 



Applying charge conjugation to an incident negative-energy particle produces 
an incident positive-energy particle, for which we know the solution of equa- 
tion ( p.55| ). Applying charge conjugation again yields the desired negative- 
energy solution. Under charge conjugation the magnitudes of the coeffi- 
cients of the scattered waves do not change, so the negative-energy solutions 
have the same reflection and transmission probabilities and as the 
positive-energy solutions. 

For the quartic Higgs potential (0), the scattering potential can 
be written in the standard form (^), so we can obtain the reflection and 
transmission probabilities Rf+ and T/-+ from equations (|104|) and (|105|) . For 
reference, the coefficients of the potential f//+ are 



^0 





Mf{l-M 



(175) 
(176) 
(177) 



where the mass ratio Aij is defined to be {^vq) / m . The scattering prob- 
abilities so obtained agree with those of Ayala et al p[. Our discussion of 



^We do not simply take w to be negative, because the analysis of asymptotic behavior 
in section 3.2 relies on the positive energy of the solution. 
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Figure 10: Transmission probability near threshold for the /+ mode, plotted 
as a function of A^|- The energy E"^ is taken to be (/if o)^ + 0.01 m^. 



fermion scattering differs from theirs primarily in that we have explicitly 
considered the polarization of the fermion. 



According to the results of section the transmission probability as a 
function of E"^ /m? has the generic form described in that section unless the 
scattering potential is close to a potential which gives T = 1, in which case 
the onset of transmission is more abrupt. Just such an increasing abruptness 
of onset occurs in the /+ mode as the mass ratio M.'j goes to zero. The effect 
can be seen in figure 0, in which we plot the transmission probability at 
a fixed distance above the threshold of total reflection as a function of the 
mass ratio Ad'j. Although the same effect occurs within the mode, it is 
of less practical importance in that case, because the gauge boson masses 
My/ and Mz are of the same order of magnitude as the Higgs mass m, while 
the fermion masses are generally much smaller. 



5 Conclusion 

We have resolved the difficulties associated with the change of particle 
content across a domain wall at a symmetry-breaking phase transition and 
obtained scattering solutions for scalar and gauge bosons. The relationship 
between the modes 02 and Bx provides a precise statement of the intuition 
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that in a spontaneously broken gauge theory each apparent Goldstone boson 
turns into the longitudinal polarization of a massive gauge boson. 

Our results should prove fundamental to many different calculations in- 
volving the electroweak phase transition in the standard model. Although 
the scattering probabilities for asymptotic particles interacting with a do- 
main wall are interesting for some purposes, we expect that the scalar internal 
modes, which are capable of describing the gradual change in particle content 
across a domain wall, will prove a more useful tool in actual calculations. 

Note added: We thank M. Voloshin for bringing to our attention earlier 
related work on this subject, in which the spontaneous breaking of a dis- 
crete symmetry was studied using a quartic potential. The single bosonic 
mode present in that situation obeys the same equations as our h mode. 
Polyakov 0] found the masses of the bound states of this mode. Voloshin 
obtained the scattering equation and hypergeometric solutions for this mode, 
including the bound-state wavef unctions, and observed that the transmission 
probability is identically one. In the latter work the scattering of fermions 
axially coupled to the scalar field was also studied. 
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